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Abstract 

SUSY23 is a Monte-Carlo package for generating supersymmetric (SUSY) processes at e+e col- 
liders. Twenty-three types of SUSY processes with 2 or 3 final state particles at tree level are included 
in version 2.0. SUSY23 addresses event simulation requirements at e+e~ colliders such as LEP. Matrix 
elements are generated by GRACE with the helicity amplitude method for processes involving massive 
fermions. The phase space integration of the matrix element gives the total and differential cross sec- 
tions, then unweighted events are generated. Sparticle widths and decay branching ratios are calculated. 
Each final state particle may then decay according to these probabilities. Spin correlations are taken into 
account in the decays of sparticles. Corrections of initial state radiation (ISR) are implemented in two 
ways, one is based on the electron structure function formalism and the second uses the parton shower 
algorithm called QEDPS. Parton shower and hadronization of the final quarks are performed through an 
interface to JETSET. 
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PROGRAM SUMMARY 

Title of program: SUSY23 

Program obtainable from: CPC Program Library, 
Queen's University of Belfast, N. Ireland (see ap- 
plication form in this issue) and from ftp.kek.jp in 
directory kek/minami/susy23. 
Computer for which the program is designed and 
others on which it is operable: HP9000 and most 
of the UNIX platforms with a FORTRAN?? com- 
piler 

Computer. HP9000 ; 

Installation: High Energy Accelerator Research 
Organization(KEK), Tsukuba, Ibaraki, Japan 
Operating system: UNIX ; 
Programming language used: FORTRAN?? 
High speed storage required: 23 Mbyte ; Card im- 
age code: ASCII 

Key words: supersymmetry, e+e~ colliders, 
ISR, QEDPS, event generator, spin correlation, 
hadronization. 



Nature of physical problem 

Study of supersymmetric particle search at LEP2. 
Method of solution 

The automatic amplitude generator GRACE is 
used to get the necessary helicity amplitudes for 
twenty-three sparticle production processes. The 
specific corrections such as the initial state radia- 
tion or hadronization are implemented in the pro- 
gram. Sparticle widths and decay branching ra- 
tios for some modes are calculated. Each event 
of final state is then generated according to these 
probabilities. 
Typical running time 

The running time depends on the number of di- 
agrams of the selected process, on the required 
cross-section accuracy and on the applied cuts. 
For instance, on a IIP-?55/99, the process e+e~ — > 
e+e/jXi takes 20 minutes to reach a 0.5% accuracy 
on the total cross section. 



1. Introduction 

The supersymmetric (SUSY) standard model is the most promising extension of the 
standard model (SM) as it could naturally give a solution to the gauge hierarchy problem. In 
SUSY model, thanks to the symmetry between bosons and fermions, the quadratic divergences 
are genuinely cancelled out. The price to pay is the existence of a large number of new particles, 
SUSY particles (sparticles), yet to be discovered. The most impressive evidence in favor of 
SUSY may be the unification of gauge couplings in SUSY Grand Unified Theories (GUTs) 
The Lightest SUSY Particle (LSP) is often proposed as a candidate to the missing mass in the 
Universe for the cold component of the so-called dark matter [^. Many peoples consider the 
"SUSY world" as a plausible scenario for future particle physics prompting the need of event 
generators. 

Although the gauge hierarchy issue set a mass upper bound to SUSY particles of O(lTeV), 
lighter sparticles are seen as giving more natural solution. So there is a possibility that the 
sparticles could be discovered at LEP2 in electron-positron collision. While the hadron colliders 
such as FNAL Tevatron and DESY HERA have already covered larger scattering energies, LEP2 
could be the first machine to discover "light" sparticles (the chargino, the slepton or the stop) 
as the background conditions are far better in electron-positron than in hadron colliders. 

SUSY23 is an event generator suitable for detector simulation and data analysis dedicated 
to SUSY processes in the LEP2 energy range. 

It is based on the GRACE ^ system which generates automatically the matrix element in 
terms of helicity amplitudes (supplied by the CHANEL |^ library) for any processes once the initial 
and the final states have been specified. In addition, a kinematics library has been developed 
for each process topology (weeding out multidimensional singularities) for a better convergence 
of the Monte-Carlo integration over the phase space. 

The SUSY23 package is actually a collection of 23 types of SUSY processes at e+e^ colliders, 
consisting of 18 final 2-body and 5 final 3-body processes, presented in a coherent and uniform 
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environment. These processes are 
e"'"e~ — > 



XiXi 


Xlx'2 


xlxl 


^2^1 










^1^2 




^R^R 




^R^L 




f^RpR 


hh 


Ul 


and 












e^e^Xi 




e^T>eXl 


^eelxi 


IXiXi 





Once a process has been selected, the total and differential cross sections are computed 
with the Monte-Carlo integration package BASES 0. Then SPRING Q, a general purpose event 
generator, provides unweighted events. Some physics results have already been presented in 
ref.i. 

The matrix element generated by GRACE corresponds to the genuine tree level process in 
the massive case. Initial state radiation corrections, unstable sparticle decays and hadroniza- 
tion must then be introduced to produce realistic event generation. All of them have been 
implemented in SUSY23. 

For the initial state radiation (ISR) two techniques are provided in the program. |^ The first 
one uses the well-known analytic form of structure function and the second is based on 
QEDPS a radiative correction generator producing an indefinite numbers of photons according 
to the parton shower algorithm in the leading-logarithmic (LL) approximation. Originally this 
algorithm has been developed to simulate QCD parton shower. One important point here is 
that QEDPS reproduces the radiative photon transverse momentum distributions. 

Sparticle widths and decay branching ratios are calculated (see Tables 3~8 in Appendix C). 
Each final state is then generated according to the branching ratios. In the present version, 
2-body and 3-body direct decays as well as some possible cascade decays are included. Note 
that helicity informations of matrix elements of the direct decays of inos ^ are taken into account 
in evaluating their decays. 

We assume that the hadronization of partons can be separated from the hard interaction 
studied here. Under this assumption, the calculation of cross sections is exact in SUSY23. Final 
state parton hadronization is performed, in SUSY23, through the mechanism implemented in 
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This paper is organized as follows. MSSM basics is given in section 2. The structure of 
the program is discussed in section 3. In section 4, all details about running the program are 
presented. A summary of the paper can be found in section 5. Four appendices describe the 
parameters and options which can be changed by the user, the list of all processes and the 
program installation procedure. 



2. Theoretical framework of MSSM 
2.1 Particle content 

^ISR is not available for the hard photon process, i.e., e~e~^ — > 7XiXi- 

^We use ino as a generic name to represent neutral (neutralinos) and charged (charginos) mixtures of 
gauginos and higgsinos. 
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The Minimal SUSY Standard Model (MSSM) includes the minimal particle content. That 
is, there should be at least one new particle (sparticle) for each known particle in the SM and 
the two Higgs doublets. The additional Higgs doublet must be included in order to give mass 
to the up- and down- type quarks and to allow for the chiral anomaly cancellation 

For quarks and charged leptons there exist two scalar partners per species, while for neutrinos 
there exists one scalar partner per species. These scalars are called squarks and sleptons, or, 
generically, sfermions. Before SU(2)xU(l) breaking, the left- and right-sfermions do not mix 
since they have different SU(2)xU(l) quantum numbers. After the breaking, however, they can 
mix with each other. In this case the mass eigenstates are parametrized by a mixing angle 9f ; 

- ( fLCOsOf + fiismef \ 
f2 J \ -Jl sin Of + Jr cos Of j 

Actually this mixing effect is substantial only for the third-generation sfermions, especially for 
the stops (superpartners of the top quark) owing to the large top mass ||l2| . 

The new fermions are either the superpartner of spin-1 gauge bosons (gauginos) or that of 
spin-0 bosons (higgsinos). They mix each other when SU(2)xU(l) symmetry is broken. The 
mass eigenstates (inos) are usually mixtures of gaugino and higgsino states. They are called 
neutralinos Xfc = 1 ~ 4) and charginos xf = 1)2), respectively, according to their electric 
charges. Gluinos g are free from mixing since the color SU(3) is not broken. As usual, we 
consider the lightest neutralino Xi ^is the lightest sparticle (LSP). 

MSSM contains two Higgs doublets and five physical Higgs bosons are left after SU(2)xU(l) 
breaking [|l3|. They are two CP even {h^ , H^), one CP odd {A'^) neutral scalars and two charged 
scalars (H^). 

The list of sparticles in the model is shown in Appendix A, where the abbreviated particle 
names used in the control card and the KF-codes assigned for each sparticles in SUSY23 are 
also shown. 

2.2 Basic parameters of MSSM 

The standard model has 18 fundamental parameters to be determined by experiments. The 
MSSM has a somewhat larger number of parameters. They are classified as (i) gauge couplings, 
(ii) superpotential parameters and (iii) soft-breaking parameters. In the following, we adopt the 
notation of the SUSY parameters described in ref . . 

2.2.1 gauge couplings 

The three gauge coupling parameters corresponding to SU(3), SU(2) and U(l) gauge groups 
are the same as in the SM. All gauge interactions are governed by these couplings. They 
determine the fermion-sfermion-gaugino interactions and four-point scalar interactions as well 
as ordinary fermion-fermion-gauge-boson interactions. 

2.2.2 superpotential parameters and i?-parity 

In the MSSM, ordinary Yukawa interactions are generalized using the superpotential W^cp), 
where (j) denotes arbitrary chiral superfield Q. Renormalizability restricts the functional form 
of the superpotential to 

W{(f)) = mij(f>i4)j + Xijk(t>i(t>j(l)k- (2) 

The parameters rriij and Xijk are further constrained by the gauge symmetry and some discrete 
symmetries. 
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Later, we use the notation /i which is a coefficient of the term quadratic in Higgs superfields 
and contributes to the mass terms of higgsinos and Higgs bosons [|T^ . 

A well-known discrete, multiplicative symmetry is the i?-parity defined by 

R=(^_fiB^L)+2S^ (3) 

where B, L and S stand for the baryon number, the lepton number and the spin, respectively. 
This formula implies that all ordinary SM particles have even i?-parity, whereas the correspond- 
ing superpartners have odd i?-parity. Usually we impose the B — L conservation on the MSSM 
and then the MSSM possesses the i?-parity invariance. In SUSY23 only i?-parity conserving 
processes are considered. 

2.2.3 soft-breaking parameters 

SUSY is not an exact symmetry of nature since our world is not manifestly supersymmetric. 
In the MSSM the SUSY breaking is induced by the soft-SUSY breaking terms, which do not 
introduce quadratic divergences. Hence the solution of the naturalness problem remains intact. 
There are four types of soft breaking terms ; (1) gaugino masses Mi (i = 1 ~ 3), (2) masses 
for the sfermions fhf, (3) trilinear term Af and (4) three Higgs mass terms. The three mass 
parameters in the last item can be re-expressed in terms of two Higgs vacuum expectation 
values, vi and V2, and one physical Higgs mass. Here vi and V2 respectively denote the vacuum 
expectation values of the Higgs field coupled to d-type and u-type quarks. vf + V2 = (246 GeV)^ 
is determined from the experimentally measured VK-boson mass, while the ratio 

tan/?=— (4) 

is a free parameter of the model. 

A generally accepted assumption is that all the three gaugino mass parameters Mj are equal 
at some grand unification scale Mx- Then the gaugino mass parameters can be expressed in 
terms of one of them, for instance, M2. The other two gaugino mass parameters are given by 

5 

Ml = -tan^ 6*^/^2 (5) 
3 

Ms = — sin^^iyMa, (6) 
a 

where a and 03 (= denote the QED and the QCD coupling constants, respectively. 

The trilinear terms Af appear in off-diagonal elements of the mass matrix for sfermions (/l, 
/ij). Consequently, it has direct relation to the mass eigenvalues {mj , fnj^ and the mixing 
angle Of. In our formulation, vl/'s (/ = t,b,T) are determined by those mixing parameters (see 
Eqs.(19), (20) and (21)). 

2.3 MSSM parameters in SUSY23 

In Appendix B, the input MSSM parameters in SUSY23 are listed. They are 

i) gaugino parameters, (tan/5, M2, /u), 

ii) masses of scalar partners of charged leptons in the 1st and 2nd generations, (^T^~^, ), 
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iii) masses of 1st and 2nd generation squarks, {jn-^, "^uh' "^5^)' ^"^ci,' '^c'h "^sij) 

iv) masses and mixing angles of charged sfermions in the 3rd generation, ijnj^^ mj^, 9t), (^^^) 

0b), (m~ , Or) and m~ . 

From the above inputs, the following physical parameters are calculated in SUSY23, 

i) masses of x° (i = 1 ~ 4), the mixing matrix On and the sign factor rfi of neutralinos, 

X^L = Vi{ON).j^^L, (7) 
X^R = r,*iONhX^R, (8) 

with ^-0, = (5^,1^0, 

ii) masses of = 1 ^ 2), the mixing angles cos^l, cos(f>R, sin^^, sin^i^, sign factor cr 
of charginos. 



XiL = W£ COS (j)L + H~^sm(pL, (9) 

X2L = -W^ sin (pL + H^L cos (l)L, (10) 

XiR = cos (t>R + H^^ sin cl)R, (11) 

X2R = eR{-Wj^ sin (pR + H^L cos (Pr), (12) 



iii) the gluino mass m~ (= M3) using Eq.(6), 

iv) m~ and m~ are calculated by the following relations: 



m~ = — mf, + m'^Y cos2P, (13) 

m~ = — + cos 2/3. (14) 



v) and are calculated by the following relations: 



vi) m~ is given by: 



^ ^i^-^u + md-mi^ cos 2/3, (15) 
rn~^ = m^^ — + — cos 2/3, (16) 



m~ = cos^ OrTn-^ + sin^ ^r'Ti- — + cos 2/3, (17) 



vii) mj_^ is given by: 



cos 6tm~ + sin 6tm~ — = cos 9bm~ + sin 6bm~ — rriu + rriuj cos 2/3, (18) 
ti 12 01 02 
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viii) the trilinear term Af using formulae as follows 

2 2 

1 mA — mt 

At = — sinlOt^ ^-^cot/3 (19) 

2 mt 

2 2 

1 mt — mt 

Ab = — sm29b^ ^-/itan/3, (20) 

2 mb 

2 2 

1 mA — mA- 

Ar = — sin20^^2 IL_^tan/3, (21) 

2 mr 

Here we briefly discuss the relation between our input parameters and those adopted in the 



susygen program |15]. As for the gaugino parameters, the sign convention of tan/3 and /x is 
different. We use the positive definite ^ and tan/? can have either positive or negative values. 
In susygen, on the other hand, positive definite tan/? and positive or negative /i are adopted. 
As there exits a symmetry for the sign of tan (3 and /.i in the MSSM, the negative ^ in susygen 
corresponds to the negative tan (3 in SUSY23. The definition of gaugino mass M2 is the same 
in both generators. As for the left-right mixing of sfermions in the 3rd generation, we take 
observable masses and mixing angles, (mj^, Of), as inputs. In susygen, on the other hand, 
more basic model parameters {mj , mj , Af) are adopted. The relation between them can be 

J— 1 

found from Eqs.(|l9|) ^ ( [2l| ) and the following formula. 



2 1 

m~ = — 

!} 2 

/2 



m~ +m~ T (("i~ - m~ )2 + (2a/m/)^) 

JL J Ft ^ JL J Ft ' 



1/2' 



(22) 



where at = —{At + cot f3), Ub = —{Ab + tan/?) and a-r = -{Aj- + tan/3). We did not impose 
low energy mass relations between masses of sfermions determined by the renormalization group 
equations in some SUSY GUTs because such mass relations sensitively depend on details of a 
model and boundary conditions for model parameters at the unification scale. In susygen, 
however, there is an option for the calculation in the framework of the minimal supergravity 



model [14|. In this model, masses of all sfermions in the 1st and 2nd generations are determined 
by an input parameter, m^o (= ttiq), which is the universal soft scalar mass at the unification 
scale. For convenience of the users, we present an example of relations between observable 
masses and the parameter moo for the case of sin^ 9w = 0.230, 

ml^ = m^ + 7.26M^ + m2 + (i-^sin2%)m|cos2/3, (23) 

= ml, + 6.83MI + ml + ^sm^9wm% cos 2P, (24) 

m~ = + 6.78M^ + - - sin^ ewml cos 2/3, (25) 

"fl 3 

= + 0.530M^ - (^ -sin2 6'H/)m|cos2/3, (26) 

= ml, + 0.1515M^+ml-sm^ewmlcos2f3, (27) 
where M^o = I.22M2 denotes the universal soft gaugino mass at the unification scale. 



3. Features of the program 



7 



This section covers some important features of the program, the introduction of the initial 
state radiation, available decay modes of unstable sparticles and the spin correlation in the ino 
decays. 

3.1 Initial state radiation 

In the first approach, the simple electron structure function is used. The electron structure 
function at O^a^) Q which is to be convoluted with the cross section for a primary process is 
given by 



D{x,s) 



U28 8 

4' ' 32 
(2a/^)(ln(s/m2)-l), 



1+3/3 + 



/3, 



1 



2-1 



\ _|_ 3j;2 

4(1 + x) ln(l - x) H In X + (5 + x) 

1 — X 



(28) 
(29) 



where s is the square of the total energy of the system and x is the momentum fraction of the 
electron. Compared with the exact 0{a) calculation for s-channel annihilation, the corrected 
cross section obtained by this function does not contain the so-called /C-factor: 




1.006480- 



(30) 



This factor is not included since it is sensible only when the s-channel e^e annihilation takes 
place. 



The basic assumption of the QED Parton Shower method, QEDPS [10|, is primarily based 
on the fact that D{x,Q^) obeys the Altarelli-Parisi equation, which can be expressed by the 
integral equation in the leading-logarithmic(LL) approximation 



D{x,Q^)=U{Q\Q^,)Dix,Ql) + ^ 

27r Jq2 



■Yi{Q\K^) 



^Piy)D{x/y,K'') 



(31) 



where the small quantity e is a cutoff related with the infrared singularity ||16|] and it will be 
defined later. In this equation P{x) is the split function noted P+{x) when regularized at x = 1. 

is the initial value of and a free parameter (of order rrig). For simplicity the fine structure 
constant a is assumed not running with Q^. The Sudakov factor 11 is given by: 



^(g^Q'^) = exp 



a 



dxP{x) 



(32) 



and denotes the probability that an electron evolves from Q'"^ to without emitting hard 
photon. 



The integral equation (31) can be solved by iteration. It is clear that the emission of n 
photons corresponds to n iterations. Hence it is possible to regard the process as a stochastic 



mechanism suggesting the shower algorithm in Ref. [10|. 

Once an exclusive process is fixed by the algorithm, each branching of a photon is a real 
process, that is, an electron with x, decays as: e~(x, —K"^) — > e^(xy, — K'^) +7(x(l — y), Qq). 
Here Q'q is a cut-off to avoid the infrared divergence and e is given by e = Ql/K'^. The 
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momentum conservation at the branching gives: —K^ = —K' /y + Qo/(l ~ u) ~\~ '^/{y{^ ~ u)) 
which in turn determines the photon transverse momentum relative to the parent, k^^, from 
y, , K''^ . This technique gives the distribution as weh as the shape of the rr-distribution. 

The above algorithm concerns the case where either e~ or e+ radiates photons when the 
axial gauge vector is chosen along the momentum of the other electron, namely e"*" or e~ . In the 
program, however, we use the double cascade scheme to ensure the symmetry of the radiation 
between e"*" and e~ |jl^. These two are mathematically equivalent in the LL approximation. 

3.2 Sparticle decays 
3.2.1 Decay modes 

The decay modes shown below are implemented in SUSY23. They are also presented in 
Tables 3~ 9 in Appendix C. The decay modes (1)~(8) are applied iteratively as long as the 
mass relationship between the sparticle and their decay products allows. The spin correlation 
is considered as is discussed in the next subsection. 

(1) lighter chargino (Table 3) 

xr ^ //'X? 



(cx?) 



(2) second lightest neutralino (Table 4) 



X2 ^ //X° 



1 



(3) left-handed selectron and smuon (Table 5) 



h ^ ^X? 

- i^Xl^i^ iffxi) 



(4) right-handed selectron and smuon (Table 6) 



^X? 

ix^^e {fM) 



(5) lighter and heavier staus (Table 7) 

n,2 TXi 

- l^rXl^'yr {fJ'xD 
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(6) sneutrinos (Table 8) 

(7) lighter stop (Table 9) 

(8) lighter sbottom (Table 9) 

6i 

We assume that squarks in the first and second generations and the gluinos are sufficiently heavy 
to forbid the charginos, neutralinos, sleptons, stop and sbottom decays in these channels. These 
are natural assumptions as long as one deal with production processes of the sparticles with 
mass ~100GeV in the LEP2 energy regions. 

In SUSY23 system, the masses are inspected whether the decay is possible or not. ^ When 
the channel is open, the decay width is computed. 

3.2.2 Spin correlations 

The helicity informations of matrix elements in the fermionic 2-body and 3-body decays of 
the lighter chargino and the second lightest neutralino X2 ^-^^ used in the event generation of 
their decay products. In other words kinematical distribution of the decay products is sensitive 
to spin correlations with the mother ino. All decay matrix elements are computed exactly with 
GRACE Q in tree approximation. 

In the event generation step, the helicity state of inos is selected by the Monte-Carlo method 
based on the relative cross sections of all possible helicity combinations. Four-momentum vectors 
of daughter-particle are determined by the decay matrix element of the ino with specified helicity. 

4. Structure of the program 

The SUSY23 provides event generation for twenty-three SUSY processes as listed in Table 2 
in Appendix C, at e^e~ colliders with or without radiative corrections. As is shown in Appendix 
B, there are many options covering theoretical and experimental requirements. Since all program 
components are distributed as source code, users can select all options by editing the relevant 
subprograms directly. However, an interface program susy23 is prepared to lighten the user's 
burden. It selects and/or corrects the program code which are affected by the various options 

^ Users must be careful in taking the input mass parameters so that the lightest neutralino Xi (LSP) 
should be lighter than charginos and sfermions. 

* Please do not confuse the executable module named susy23 with the SUSY23 system itself. 



(//X?) 
(//'X?) 



CX? 
CX^ 

&xr 



c (//X?) 
>h (Z/'x?) 



&X? 
6X^ 



h iffxl) 
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and it creates a "Makefile" according to the user requirements. This procedure is called the 
source generation phase. 

In the integration step, the matrix element of a selected process is integrated over the phase 
space by the subprogram BASES, which gives the total and differential cross sections and the 
probability distribution used in the event generation phase Q. There, the subprogram SPRING 
samples a point in the phase space and test if it can be accepted as a new event according to 
its probability. When an event is accepted, the program control returns to the main program, 
where sparticle decays are performed. 

There are, therefore, three steps in the generator SUSY23, the first is the source generation, 
the second is the integration and the third is the event generation. In addition to the user 
interface program susy23, the following programs are available in the SUSY23 system: 

i) Twenty-three function programs FUNCs for final 3-body or 2-body SUSY processes, each 
of which calculates the numerical value of the differential cross section for each process. 

ii) The kinematics subprogram KINMOQ is called when initial state radiation (ISR) based on 
QEDPS, the QED parton shower model, is requested. Although KINEMO covers both the 
case where no ISR is produced and where at most 2 photons can be generated at zero 
angle following the electron structure function method. 

iii) An extended version of the CHANEL ^ library to include MSSM couplings. The CHANEL 
library calculates the numerical values of Feynman diagrams based on the helicity ampli- 
tudes method. 

iv) The numerical integration and event generation program package BASES/SPRING v5 . 1 [Q]. 

v) The main programs, MAINBS and MAINSP, and all program components for the integration 
and event generation steps. 

vi) The subprograms for branching ratios and sparticle decays. 

v) The interface programs to JETSET (GR2LND, SP2LND and GRC2SH). 

The function of these program components and relationship among them are presented in the 
next three subsections. 

4.1 Source generation step 

The user interface program susy23 reads the parameters from the control data file, which 
contains process selection, type of radiative corrections, SUSY parameters, etc. The complete 
list of parameters is found in Appendix B. Then susy23 generates the following program code: 

i) Three initiahzation subprograms USRPRM, MODMAS KINIT, USERSP and 

ii) A "Makefile". 

4.2 Integration step 

Before starting the numerical integration, the main program MAINBS invokes an initialization 
subprogram USERIN, in which the following subprograms are called in this order: 
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USRPRM : To define the set of input SUSY parameters and some optional 
parameters. 

SETMAS : To set default values to particle masses and decay widths. 
MODMAS : To alter the default values of all parameters defined in SETMAS. 
AMP ARM : To set the coupling constants and other parameters. 
KINIT : To set the parameters for the integration, kinematics, cuts and to 
initialize the histogram package. 

The subprograms SETMAS and AMPARM are generated by the GRACE system. By default ino 
and charged Higgs masses, as well as ino mixing angles (see Sec. 2. 3) and all branching ratios, 
are calculated in the framework of the MSSM. Users can selectively set the values of these 
parameters and some branching ratios. Only naive consistency checks are made, such as the 
normalization check of the branching ratios and the mass ordering check, so all modifications 
on these parameters or branching ratios in the subprograms USRPRM, MODMAS and KINIT are on 
the users' own responsibility. 

The integration program BASES calculates the scattering cross section by sampling the func- 
tion FUNC on the allowed phase space segmented by a self adapted grid where finer cells are 
clustered on the high gradient zones. This is an iterative procedure running until either the 
maximum number of allowed iteration is reached or the required accuracy is obtained. In the 
function program FUNC, the kinematics subroutine KINEMO or KINMOQ maps the integral vari- 
ables into the four-momentum of the final state particles. KINEMO is called for reactions with 
no radiative corrections or those involving the initial state radiation treated with the structure 
function techniques. KINMOQ is called for processes in which radiative corrections are computed 
with the QED parton shower method. The subprograms AMPTBL and AMPSUM are further called 
for calculating the helicity amplitudes and the squared of their sum. 

It is recommended to look at the integration result carefully, especially over the 
convergency behaviors both for the grid optimization and integration steps. When 
the accuracy of each iteration fluctuates from iteration to iteration or when it jumps up suddenly 
to a large value compared to the other iterations, the resultant estimate of the integral may not 
be reliable. There are two possible origins of this behavior; too few sampling points or/and an 
unsuitable choice of the kinematical variables. 

After the numerical integration by BASES, the subprograms BSINFO and BHPLOT are called 
to print the result of integration and the histograms, respectively. Before terminating the in- 
tegration procedure the probability distribution of the integrand can be saved in a file named 
bases. data by invoking BSWRIT, which is later used for the event generation by SPRING. 

4.3 Event generation step 

The main program, MAINSP, handles the event generation step. The subprogram BSREAD 
reads the data file bases .data to restore the probability distribution and then the subprogram 
USERIN is called. Each call to SPRING generates one event by sampling a point in the phase 
volume. It first calculates the differential cross section at that point using the same function 
FUNC in the integration step and returns the weight of this sampling point. A weight one event is 
finally produced using usual unweighting technique. Subsequently, SUDCAY handles the sparticle 
decay. By calling the subprogram SP2LND, the event information including the color connection 
to be used by JETSET is stored in the labelled common LUJETS. 

Then hadronization of quarks and gluons can be performed by calling LUEXEC. At the end of 
the event generation, the routine SPINFO and SHPLOT are invoked successively for printing event 
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generation information and histograms. 



5. How to run the program 

The user should first prepare the control data to define the process, the option flags, the 
MSSM parameters and the experimental cuts. The user interface program susy23 takes this 
control data as an input. For example, let's call the following sequence of lines control. data. 

process = SWlswl 

energy = 200 . OdO 

type = tree 

hadron = no 

tanbe = -2 . OdO 

xiii2 = 50 . OdO 

xmu = 150 . OdO 

end 

The first line specifies the process to be calculated and the second is the center of mass energy 
in GcV unit. The others are options, whose meanings are given in Appendix B. Then the user 
may type: 

% susy23 < control. data 

If the message "syntax error" is returned, the user must carefully examine the contents of the 
control data. No file is generated in this case. When the control data is accepted the following 
message appears: 

Process is "SWlswl" 
Energy is "200. OdO" 
HADRNZ <no> 
bye-bye 

directory name is SWlswl 

absolute directory name is /home/susy23/prc/SWlswl 

cd /home/susy23/prc/SWlswl 

make integ 

integ 

make spring 
spring 



According to the parameters given in the control data, the files, i.e. usrprm.f, modmas.f, 
kinit.f, usersp.f and Makefile, are generated in a specified subdirectory ( SWlswl in this 
case). 

According to the last five lines in the message, users can proceed with the calculations as follows: 
i) Change directory by typing: 
% cd /home/susy23/prc/SWlswl 



13 



ii) Create an executable integ for the integration by typing: 
% make integ 

iii) Numerical integration is actually performed by typing: 
y. integ 

The results of integration step arc displayed on the console as well as written in an out- 
put file bases . result. The total cross section in pb and the estimated statistical error 
are shown on the last line, under Cumulative Result, in the table of the Convergence 
Behavior for the Integration step. The differential cross sections are also printed as 
a function of the energy, scattering angle of each particle and invariant masses of any two 
final particles. The probability distribution of the integrand is written in a file bases . data 
which will be used in the event generation step by spring. 

iv) Before running the event generation, users may edit mainsp . f to set additional parameters 
if needed and call user's own analysis routines. 

The following is the structure of the generated mainsp. f, where four- momentum of all 
particles arc stored in the common/lujets/ in the JETSET format when subprogram sp21nd 
is called in the event-loop: 

Program mainsp 
implicit real*8(a-h,o-z) 
external func 



real*4 p,v 

common/lujets/n,k(4000,5) ,p(4000,5) ,v(4000,5) 



mxtry = 50 
mxevnt = 1000 
do 100 nevnt = 1 , mxevnt 

call spring ( func, mxtry ) 



* 

call sp21nd 

* 

* 

* ============================================= 

* ( user_analysis based on the common lujets ) 
Mc ============================================= 

* 

100 continue 



stop 
end 
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v) Create an executable spring for event generation by typing: 
% make spring 

vi) Start the event generation by typing: 
°/„ spring 

Information on the event generation wih be written in the spring. result file. Users 
should pay attention to the histograms generated in this step. The distributions of the 
generated events are superimposed with the character "0" on the histograms generated in 
the integration step. These two distributions should be consistent with each other within 
the statistical error of the generation. For the detail of the output files of BASES and 
SPRING, users can consult the Ref.Q. 



6. Summary 

The SUSY23 system calculates the effective cross section and generates events for one of the 
twenty-three SUSY processes at e^e~ colliders listed in Appendix C. It is dedicated to the LEP2 
studies. The numbers 2 and 3 in SUSY23 refer to the fact that processes with 2 and also 3 final 
particles are fully computed at tree level in the massive case. These particles are then allowed 
to decay producing multi-particle final states. Sparticle widths and decay branching ratios in 
accordance with this energy range are calculated. In particular, the helicity information for 
direct decays of spinor sparticles (the lighter chargino and the second lightest neutralino) is 
handled properly. The produced quarks can be hadronized according to JETSET. Processes with 
initial radiations can be generated in terms of the electron structure function or the QED parton 
shower method. 
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Appendix A. List of sparticles in the MSSM 

The list of sparticles in the MSSM is shown in Table 1, which includes abbreviated names 
used in the control cards and the KF-codes asigned for each particles. As is seen, the keywords 
in SUSY23 are case-sensitive. 



particle 


abbrev. 


KF code 


particle 


abbrev. 


KF code 


xt 


SWl 


75 




SUL 


-42 


xl 


swl 


-75 


UR 


sur 


48 


xt 


SW2 


76 


^R 


SUR 


-48 


X2 


sw2 


-76 


dL 


sdl 


41 


X? 


szl 


71 


dL 


SDL 


-41 


x'2 


sz2 


72 


dR 


sdr 


47 


xt 


sz3 


73 


dn 


SDR 


-47 


xl 


sz4 


74 


CL 


scl 


44 


9 


sgl 


70 




SCL 


-44 


eZ 


sel 


51 


cr 


scr 


50 


et 


SEL 


-51 


^R 


SCR 


-50 




ser 


57 


SL 


ssl 


43 




SER 


-57 




SSL 


-43 




sne 


52 


■V-!. 


ssr 


19 




SNE 


-52 


^R 


SSR 


-49 


AZ 


sml 


53 


ii 


stl 


46 




SML 


-53 


h 


STl 


-46 




smr 


58 


h 


st2 


62 


I^R 


SMR 


-58 


h 


ST2 


-62 




snm 


54 


bi 


sbl 


45 




SNM 


-54 


bi 


SBl 


-45 




sal 


55 




sb2 


61 


ft 


SAl 


-55 


h 


SB2 


-61 




sa2 


59 




shl 


25 


ft 


SA2 


-59 




sh2 


35 




sna 


56 


AO 


sh3 


36 




SNA 


-56 


H+ 


sh4 


37 


UL 


sul 


42 


H- 


SH4 


-37 



Table 1 Sparticles in SUSY23. 



Appendix B. Parameters in control data 

In the table below, the default values are underlined and the relation between commands and 
variable/array in Fortran sources is also described. Variable names are written in bold letters 
and filenames in italic. All masses and energy parameters are given in GeV. 
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i) Process selection. 



Process = SWlswl 

abbreviation of process name 

This specifies the subdirectory name, where subroutines for the process are stored. Table 
1 in Appendix C shows the abbreviation of process names. 

ii) Center of mass energy. 

energy = 200 .DO 

CMS energy 
w in kinit.f 

iii) Control parameters. 



helicityl = 


average, left, right 




Helicity state for the initial electron. 


helicity2 = 


average, left, right 




Helicity state for the initial positron. 



type 


= tree, sf , qedpsi (Not available for ^yYiYl) 




Type of calculation: 




Without radiation(tree), ISR with structure function(sf ) 




and ISR with QEDPS(qedpsi). 




jqedps in usrprm.f, isr in kinit.f. 




jqedps = 0, isr = without radiation(tree). 




jqedps = 0, isr = 1 for sf . 




jqedps = 1 for qedpsi. 



decay 


= yes, no 

Unstable sparticle decays : 
j dec ay = 1, or in usrprm.f. 


spincr 


= yes, no 

Include spin correlations for direct 3-body decays of gaug- 
inos : 

j spncr = 1 or in usrprm.f. 


qcdcr 


= yes, no 

Include overall QCD correction factor in squark produc- 
tion : 

j qcdcr = 1 or in usrprm.f. 


susyp 


= yes, no 

(yes) user can selectively set ino masses and mixing an- 
gles as well as charged Higgs mass in modmas.f (see next 
sections) 

(no) all parameters are calculated in the MSSM frame- 
work : 

j susyp = 1 or in usrprm.f. 



iv) Standard Model parameters. 
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v) 



amw 


= 80 . 23D0 

W mass : amw in modmas.f. 


agw 


= 2.08D0 

W width : agw in modmas.f. 


amz 


= 91.19D0 

Z mass : amz in modmas.f. 


agz 


= 2.491D0 

Z width : agz in modmas.f. 


amu 


= lOO.OD-3 

u-quark mass : amu in modmas.f. 


amd 


= lOO.OD-3 

d-qiiark uimss : amd in modmas.f. 


amc 


= 1.8645D0 

c-quark mass : amc in modmas.f. 


ams 


= 300.0D-3 

s-quark mass : ams in modmas.f. 


amt 


= 174. ODO 

top mass : amt in modmas.f. 


amb 


= 5 . ODO 

b-quark mass : amb in modmas.f. 




alphai 


= 128.07D0 

a^^ : alphai in modmas.f. 


alpha_s 


= . 12D0 

as : alphas in modmas.f. 


Minimal SUSY Standard Model parameters. 


tanbe 


= -2 . ODO 

tan/3 : tanbe in usrprm.f 


xm2 


= 50. ODO 

SU(2) gaugino mass M2 : xm2 in usrprm.f. 


xmu 


= 150. ODO 

SUSY Higgs mixing mass fi : xmu in usrprm.f. 
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Sfermion sector parameters. 



XI1IS61 


= 95.0D0 

left-handed selectron mass m — : xmsel in usrprm.f. 


xmssr 


= 90.0D0 

right-handed selectron mass m — : xmser in usrprm.f. 


xmsml 


= 95 . ODO 

left-handed smuon mass m — : xmsml in usrprm.f. 


xmsmr 


= 90. ODO 

right-handed smuon mass m~- : xmsmr in usrprm.f. 

Mr 




xmsal 


= 90. ODO 

lighter stau mass m~ : xmsal in usrprm.f. 


xmsa2 


= 95 . ODO 

heavier stau mass m~ : xmsa2 in usrprm.f. (this param- 
eter is recalculated in case of susyp = no) 


xmsna 


= 100. ODO 

tau-sneutrino mass : xmsna in usrprm.f. 


tha 


= pi/2. ODO 

mixing angle Or of staus : tha in usrprm.f. 




xmsul 


= 200. ODO 

left-handed s u-quark mass : xmsul in usrprm.f. 


xmsur 


= 200. ODO 

right-handed s u-quark mass rriy^^ : xmsur in usrprm.f. 


xmsdr 


= 200. ODO 

right-handed s d-quark mass m j : xmsdr in usrprm.f. 

an 


xmscl 


= 200. ODO 

left-handed s c-quark mass : xmscl in usrprm.f. 


xmscr 


= 200 . ODO 

right-handed s c-quark mass ^n^^ : xmscr in usrprm.f. 


xmssr 


= 200. ODO 

right-handed s s-quark mass rrig^ : xmssr in usrprm.f. 
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xmstl 


= 80.0D0 

iignter stop mass mj^ : xmstl m usrprm.j. 


xmst2 


= 250. ODO 

iKVuicr st()|") uiMSs irij^ : xmst2 iu iisrpnn.f. 


xmsbl 


= 90. ODO 

lighter sbottom mass m^^ : xmsbl in usrprm.j. 


xmsb2 


= 210. ODO 

heavier sbottom mass m^^ : xmsb2 in usrprm.f. This pa- 
rameter is calculated in case of susyp = no. 


tht 


= 1 . 5D0 

mixing angle Of of stops : tht in usrprm.f. 


thb 


= . 3D0 

mixing angle 9b of sbottoms : thb in usrprm.f. 



xaa 


= . ODO 








trilinear coupling Ar 


in modmas.f in case of susyp 


= yes. 


xau 


= 250 . ODO 








trilinear coupling At 


in modmas.f in case of susyp 


= yes. 


xad 


= 50. ODO 








trilinear coupling A^ 


in modmas.f in case of susyp 


= yes. 
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vii) Neutralino sector parameters. 

For only susyp = yes, the following parameters can be given by users. 

If susyp = no is specified, the input for the following parameters are neglected since these 

arc calculated by Eqs.(7) and (8). 



amnl 


= 30 . ODO 

neutralino mass Xi- 


aiim2 


= 70 . ODO 

neutralino mass X2- 


amnS 


= 200. ODO 

neutralino mass Xs- 


aim4 


= 400 . ODO 

neutralino mass Xi- 


on 


= l.ODO ffor i= 1 4) 

= O.ODO (for i/j, i,j = 1 ~4) 

neutralino mixing matrix Ctv in modmas.f 


sgn 


= i, l.ODO, 0. ODO ffor i = 1 ~ 4) 

neutralino sign factor Re{r]i), Im{r]i). 



viii) Chargino sector parameters. 

For only susyp = yes, the following parameters can be given by users. 

If susyp = no is specified, the input for the following parameters are neglected since these 
are calculated by Eqs.(9) ~ (12). 



amcl 


= 30 . ODO 

chargino 


mass xr- 




aiiic2 


= 200. ODO 

chargino 


mass X2 • 




cphl 


= 1 . ODO 








chargino 


mixing matrix cos 


h- 


sphl 


= . ODO 








chargino 


mixing matrix sine; 


i>L. 


cphr 


= 1 . ODO 








chargino 


mixing matrix cos 




sphr 


= . ODO 








chargino 


mixing matrix sine; 




eps 


= 1 . ODO 

chargino 


sign factor cr. 





ix) Experimental cuts. 

The numbering convention of particles follows the GRACE scheme, where the initial electron 
and positron are 1st and 2nd, respectively, and the final particles are numbered 3, 4 and 
5. In the process name of Table 2, the order of particles corresponds to this numbering 
convention. For instance, in the process, e+e^Xi ) the 3rd is , the 4th is and the 5th 
is Xi- 
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thecutS 


= 180. DO, O.ODO 






Polar angle cut for 3rd particle in degree (backward-anj 






forward-angle) . 






coscut (1 : 2 , l)=cos (thecutS) in kinit.f. 




thecut4 


= 180. DO, O.ODO 






Polar angle cut for 4th particle in degree (backward-anj 






forward-angle) . 






coscut (1 : 2 ,2)=cos (thecut4) in kinit.f. 




thecutS 


= 180. DO. O.ODO 






Polar angle cut for 5th particle in degree (backward-anj 






forward-angle) . 






coscut (1 : 2 ,3)=cos (thecutS) in kinit.f. 





Instead of giving a numerical value, the user can use the strings as below: 
amassl(i) has the mass for i-th particle and w is the CM energy. 



engcutS 


amasslO), w 




Energy cut for 3rd particle (min.,max.) 




engyct (1 : 2 , 1) in kinit.f 


engcut4 


amassl(4), w 




Energy cut for 4th particle (min.,max.) 




engyct (1 : 2 ,2) in kinit.f 


engcutS 


amassl(5), w 




Energy cut for 5th particle (min.,max.) 




engyct (1 : 2,3) in kinit.f 



x) Parameters for integration. 



itmx 


= 5,10 




Iteration numbers: itmxl, itmx2 in kinit.f 


acc 


= 0.1, 0.05 




Accuracies in %:accl, acc2 in kinit.f 


ncall 


= 10000 




Sampling points: ncall in kinit.f 



xi) Parameters for event generation. 



mxtry 



50 

Maximum trial numbers: mxtry in usersp.f 



mxevnt 



10000 

Maximum event numbers: mxevnt in usersp.f 



hadron 



yes, no 

Hadronization with JETSET|11| for event generation: 

If yes, then a statement, call luexec is activate in 

grc2sh.f. 



xii) End of description. 



end 



After the command end any command is neglected. 
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Appendix C. Process table 

The 23 processes included in SUSY23 are listed in Table 2. In the heading, 'abbrev.' and 
'dir.' stand for the abbreviated name used in the control card and directory name where the 
generated code is stored, respectively. In Tables 3 ~ 9, decay modes taken into account in the 
unstable sparticle decays in SUSY23 are listed. By default, each branching ratio is automatically 
calculated in the framework of the MSSM. The abbreviated names are used in the control card 
or kinit.f when users want to set each branching ratio (jbrnch = 1). 



process 


abbrev. 


dir. 


process 
e'^e~ — > 


abbrev. 


dir. 


xtxi 


SWlswl 


SWlswl 




SA2sal 


SA2sal 




szlsz2 


szlsz2 




SA2sa2 


SA2sa2 


x^x^ 


sz2sz2 


sz2sz2 




SNEsne 


SNEsne 


TXiXi 


aszlszl 


szlszla 




SNMsnm 


SNMsnm 


^R^R 


SERser 


SERser 




SNTsnt 


SNTsnt 




SELsel 


SELsel 


e^e]^.Xi 


Eserszl 


Eserszl 


^P^R 


SELser 


SELser 


e+e^Xi 


Eselszl 


Eselszl 


^R^L 


SERsel 


SERsel 


e^J^eXi 


Esneswl 


Esneswl 


f^Rp^R 


SMRsmr 


SMRsmr 


^e^hXl 


nSELswl 


nSELswl 




SMLsml 


SMLsml 


tJl 


STlstl 


STlstl 




SAlsal 


SAlsal 




SBlsbl 


SBlsbl 




SAlsa2 


SAlsa2 









Table 2 Available processes in SUSY23. 



process 


abbrev. 


process 


abbrev. 


Xl e-UeXi 


swlbr(l) 


Xr '^"^r 


swlbr(8) 


Xl l^'^fiXi 


swlbr(2) 


Xl T^eel 


swlbr (9) 


Xl t'^tXi 


swlbr(3) 


Xl ^^fil 


swlbr(lO) 


Xl udxi 


swlbr(4) 


Xl Vrfi 


swlbr(ll) 


Xl csxi 


swlbr (5) 


xr ^rT2 


swlbr(12) 


Xl e'Ue 


swlbr (6) 


Xl bil 


swlbrdS) 


Xl l^'^iJ, 


swlbr (7) 







Table 3 Decay processes of lighter chargino Xi in SUSY23. 



25 



process 


abbrev. 


process 


abbrev. 


Yo — > e~e^Y? 


sn2br(l) 


A2 1 


sn2br(16) 


Yo — > u~u^y^ 


sn2br (2) 


yO ^ 

A.2 1 


sn2br(17) 


Yo — > T~T~'~Y? 


sn2br (3) 


Yo — >■ e^e7 

A2 ^ L 


sn2br(18) 


Yo — VcUfY^ 
A2 '^e'^eAl 


sn2br (4) 


Yo — > e~ et 

A2 ^ L 


sn2br(19) 


Yo — > ViiUlY^ 
A 2 *^/i*^/iAl 


sii2br (5) 


Yo — > Ll^ U, T 
A 2 A L 


sn2br (20) 


Yo — * V i' 

A2 ' / ' / Ai 


sn2br (6) 


-V 2 / ' 1 L 


sn2br (21) 


Yo — > li'UY? 
A2 ^Al 


sii2br (7) 


A2 ' ' 2 


sn2br (22) 


yO ^ 


sn2br(8) 


A2 ' ' 2 


sn2br (23) 


Yo — > CCY? 


sn2br (9) 


Yo — !■ Z^oP^ 
A2 *^e^e 


sn2br(24) 


Yo — > 5Sy9 
A2 '-^'-^Al 


sn2br(10) 


Yo — !■ t'eP^ 
A2 '^e'^e 


sn2br (25) 


yO _^ 55 
Az ^^Al 


sn2br(ll) 


Yo — > ^'//^';, 
A2 fl^ fl 


sn2br (26) 


X2 ^ e+e;^ 


sn2br(12) 


X2 ^H^H 


sn2br(27) 


X2 ^ e-e+ 


sii2br(13) 


X2 ^tI>t 


sn2br(28) 




sn2br(14) 


A:2 — ''r'^^r 


sn2br(29) 


X2 /^"Ar 


sn2br(15) 







Table 4 Decay processes of second lightest neutralino xi SUSY23. 



process 


e = e 


e = n 


h^^xl 


selbr(l) 


smlbr(l) 


h^^xl 


selbr(2) 


smlbr(2) 


h J^exT 


selbrO) 


smlbr(3) 



Table 5 Decay processes of left-handed charged sleptons £l in SUSY23. 



process 


i = e 


t = li 


h^^xl 


serbr(l) 


smrbr ( 1 ) 


Ir ^ H-^2 


serbr (2) 


smrbr (2) 



Table 6 Decay processes of right-handed charged sleptons Ir in SUSY23. 
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process 


abbrev. 


Ti~ — > T~y9 

'1 ' Al 


salbr(l) 


Tt~ — >■ T~Yn 


salbr(2) 




salbrO) 


^2 ^"Xl 


sa2br(l) 


T2 T'X2 


sa2br(2) 


T2 /^r\T 


sa2br(3) 



Table 7 Decay processes of lighter and heavier staus f in SUSY23. 



process 


£ = e 


e = n 


e = T 




snebr (1) 


smnbr (1) 


snabr (1) 


h ^tX2 


snebr (2) 


smnbr (2) 


snabr (2) 


h (-xt 


snebr (3) 


smnbr (3) 


snabr (3) 



Table 8 Decay processes of sneutrinos in SUSY23. 



process 


abbrev. 


ii ^ cx? 


stlbr(l) 


ii ^ cx^2 


stlbr(2) 


h ^ bxt 


stlbr(3) 


bi ^ 6x? 


sblbr(l) 


h - bxl 


sblbr(2) 



Table 9 Decay processes of 3rd generation squarks ti and bi in SUSY23. 

Appendix D. Installation 

The source code is available by anonymous ftp from ftp.kek.jp in the directory 
kek/minaiiii/susy23. The SUSY23 system contains the complete set of Fortran sources for 23 
processes, the three libraries, i.e., BASES/SPRING, CHANEL and utilities for kinematics. Those 
source codes are written in FORTRAN77. In addition, SUSY23 provides the interface program 
to generate a few Fortran source files according to the control data specified by the user. This 
program is written in C, YACC and LEX. SUSY23 has been developed on HP-UX, but should 
run on any UNIX platform with a fortran compiler. 

The procedure of installation is as follows: 

1. Editing Makefile. 

The following macros in Makefile should be taken care of by users themselves. For example, 
SU23DIR defines the directory name where SUSY23 is installed. The values of FC and FOPT 
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define the relevant compiler name and option for your system. The other macros can be 
left as they are. 

SU23DIR = directory where SUSY23 are installed. 
PRCDIR = directory where process files are installed. 

(default is $(SU23DIR)/prc.) 
LIBDIR = directory where libraries are installed. 

(default is $(SU23DIR)/lib.) 
BINDIR = directory where an executable is installed. 

(default is $(SU23DIR)/bin.) 
MACHINE = [hpux | hiux | sgi | dec | sun] 
FC = FORTRAN compiler command name. 

FOPT = FORTRAN compiler options. 

CERNLIBS = CERNLIB including the jetset library. 

For HP-UX, -L/cern/pro/lib -ljetset74 

-Ipacklib -Ikernlib -L/lib/pal . 1/ -Im 

2. Compilation. 

By executing command make install the executable of the interface program(susy23) is 
generated at BINDIR. Furthermore six libraries, BASES/SPRING, CHANEL, kinematics utility 
library and three susy libraries, arc generated in LIBDIR. 



The sample control data files will be found in the directory sample. 



TEST RUN OUTPUT 
control data 



Process 


= Eselszl 


Energy 


= 200. OdO 


type 


= tree 


itmx 


= 7, 15 


ncall 


= 10000 


hadron 


= yes 


mxevnt ■ 


= 10000 


end 





Followings are the output files from BASES and SPRING. Only one histogram, the momentum 
distribution of the particle 3, is shown since the whole output is too lengthy to be included here. 
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— ####-##-YYYY- 
-########— YYYY- 

-## ###— YY- 

-####### YY- 



222 

2222222— 

22 22- 

-YYYY 22- 

-YYYY 22- 

-YY 22— 

-YY 222— 



* SSSS SS-UUU— UUUU— 

*— ssssssss— uuu— uuuu— 

*— ss sss— uu uu— 

*— sssssss uu uu— 

* sssssss— uu uu #######— YYY-YY 222 

*— SS SS— UU— uuu— ## ## YYYYY 22 22- 

*— SSSSSSSS UUUUUUUU— ######## YYYY 22222222- 

*— SSSSSSS UUU-UUU— ####### YY 22222222- 

* YYY 

* YYYYYYY 

* YYYYYYY 



333 * 

—3333333 * 

—3 33 * 

33 * 

3333 * 

3333 * 

33 * 

33 * 

—3 33 * 

-3333333 * 

333 * 



* matrix element generated by GRACE * 

* version 2.0 (07) — 1997 JUN. 17 * 

* Copyright — Mineuni Tateya Collaboration * 

* E-mail : susy23aminami .kek. jp * 

HcHcHc^Hc^Hc^Hc^^Hc^Hc^Hc^Hc^Hc^H^^st^H^^^^^^^^Hc^^Hc^^Hc^HcHc^HcHc^Hc^^tEHc***** ************ 

* Process : e+ 



CM Energy (GeV) : 

Mass (Width) 

W-boson : 
u-quark : 
c-quark : 



— > (3)positron (4) selectronl (5)neutralinol 
200.00000 



80. 230 ( 
.100 
1.865 
t-quark: 174.000 



2.080) 



2-boson: 
d-quark : 
s-quark: 
b-quark : 



91. 190( 

.100 

.300 
5.000 



2.491) 



alpha : 
alpha_s : 



1/128.070 
. 12000 



SUSY Parameters 



tan (beta) 

SU(2) gaugino mass(Gev) 
SUSY Higgs parameter (Gev) 



-2.00 
50.00 
150.00 



SFermions Masses (GeV) 



s-electron(L) : 


95, 


,000 


s-electron(R) : 


90, 


,000 






s-neutrino-electron : 


95, 


.000 


s-muon(L) : 


95, 


.000 


s-muon(R) : 


90, 


.000 








s-neutr ino-muon : 


95, 


.000 


s-tau(l) : 


90, 


.000 


s-tau(2) : 


95, 


.000 








s-neutr ino-tau : 


100, 


.000 


s-u-quark(L) : 


200. 


,000 


s-u-quark(R) : 


200, 


.000 


s-d-quark(L) : 


200. 


,000 


s-d-quark(R) : 


200, 


.000 


s-c-quark(L) : 


200. 


,000 


s-c-quark(R) : 


200, 


.000 


s-s-quark(L) : 


200. 


,000 


s-s-quark(R) : 


200, 


.000 


s-t-quark(l) : 


80, 


.000 


s-t-quark(2) : 


250, 


.000 


s-b-quark(l) : 


90, 


.000 


s-b-quark(2) : 


210, 


.000 



SFermions Mixing Angles (rad.) 



s-tau: 
s-top: 



.000 
L.500 



s-bottom: 



.300 
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trlinear couplings 




A_tau : 


.3000000E+03 


A_u: 


. 5225048E+02 


A_d: 


- . 1732713E+04 



Neutral Ino 



neutralino masses (Gev) 



27.89, 
159.24, 



70.53 
183.34 



* 


neutralino mixing 


matrix (bino 


bases) 




* 


.952 


.284 


104 


.036 


* 


-.207 


.864 


459 


-.018 


* 


-.122 


.280 


611 


-.731 


* 


-.187 


.308 


637 


.681 



Neutralino sign factor 
1.00 .00 I 1.00 



.00 I 1.00 



.00 I 



.00 1.00 



Chargino 

chargino masses (GeV) 
chargino mixing matrix 

chargino sign factor 



69.60, 



.85, 
1.00, 

1.00 



181.74 



-.53 
.09 



* 
* 


Calculated Total Decay width(GeV) 






* 


chargino 1 


.364E-04 


charglno2 


.521E+01 




* 


neutralino2 


. 149E-04 


neutralinoS 


.322E+00 neutralino4: 


. 662E+00 


* 


selectronR 


. 339E+00 


selectronL 


.392E+00 




* 


smuonR 


. 339E+00 


smuonL 


. 392E+00 




* 


staul 


.325E+00 


stau2 


. 367E+00 




* 


sneu_e 


.272E+00 


sneu.mu 


. 272E+00 sneu.tau : 


.351E+00 


* 


stopl 


.240E-03 


sbottoml 


.469E-01 





Experimental Cuts 
Angle Cuts 



* 


Particle 3 


-1 


000 


> 


1 


000 


* 


Particle 4 


-1 


000 


> 


1 


000 


* 


Particle 5 


-1 


000 


> 


1 


000 


* 


Energy CUT 












* 


Particle 3 




001 


> 


62 


246 


* 


Particle 4 


95 


000 


> 


120 


618 


* 


Particle 5 


27 


888 


> 


79 


382 


* 


Invariant mass CUT 












* 


Inv. Mass 4-5 


122 


888 


> 


199 


999 


* 


Inv. Mass 3-4 


95 


001 


> 


172 


112 


* 


Inv. Mass 3-5 


27 


889 


> 


105 


000 



* OPTIONS: * 

* Calculation: TREE * 

* Width : FIX * 
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Branching Ratios of Sparticles 



chargino 1 


(swl) 










(EL ne szl) 


.1774 


(MU nm szl) 


.1774 


(TA nt szl) 


.2171 


(uq DQ szl) 


.2153 


(cq Sq szl) 


.2129 






(EL sne) 


.0000 


(MU snm) 


.0000 


(TA snt) 


.0000 


(ne SEL) 


.0000 


(nm SML) 


.0000 






(nt SAD 


.0000 


(nt SA2) 


.0000 






(BQ stl) 


.0000 










neutralino 2 


(sz2) 










(EL el szl) 


.2747 


(MU mu szl) 


.2747 


(TA ta szl) 


.3367 


(NE ne szl) 


.0200 


(NM nm szl) 


.0200 


(HT nt szl) 


.0136 


(UQ uq szl) 


.0198 


(DQ dq szl) 


.0017 






(CQ cq szl) 


.0194 


(SQ sq szl) 


.0017 


(BQ bq szl) 


.0177 


(el seR) 


.0000 


(mu smR) 


.0000 


(ta sal) 


.0000 


(el seL) 


.0000 


(mu smL) 


.0000 


(ta sa2) 


.0000 


(ne sne) 


.0000 


(nm snm) 


.0000 


(na sna) 


.0000 


selectron L 


(seL) 










(el szl) 


.5581 


(el sz2) 


.1193 


(ne swl) 


.3226 


selectron R 


(seR) 










(el szl) 


.9914 


(el sz2) 


.0086 






smuon L 


(smL) 










(mu szl) 


.5581 


(mu sz2) 


.1193 


(nm swl) 


.3226 


smuon R 


(smR) 










(mu szl) 


.9914 


(mu sz2) 


.0086 






stau 1 


(sal) 










(ta szl) 


.6223 


(ta sz2) 


.1007 


(na swl) 


.2770 


stau 2 


(sa2) 










(ta szl) 


.9883 


(ta sz2) 


.0115 


(na swl) 


.0002 


snu_e 


(sne) 










(ne szl) 


.0669 


(ne sz2) 


.2900 


(el swl) 


.6431 


snu.mu 


(snm) 










(nm szl) 


.0669 


(nm sz2) 


.2900 


(mu swl) 


.6431 


snu.tau 


(sna) 










(nt szl) 


.0555 


(nt sz2) 


.2960 


(tu swl) 


.6485 


stop 1 


(stl) 










(cq szl) 


.0000 


(cq sz2) 


.0000 


(bq swl) 


1 . 0000 


sbottom 1 


(sbl) 










(bq szl) 


.1410 


(bq sz2) 


.8590 







>|c * >|c :1c >|c * >|c * >|c * >|c :1c >|c * >|c 3|c >|c * >|c :1c >|c * * >|c :1c >|c >|c * >|c >|c 3|c * 3|c 3|c :|c 3|c 3|c :|c 3^ 
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AA 


AA 
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* 


* 
* 






BASES 


Version 5 


1 




* 
* 


* 


coded by S.Kawabata KEK 


March 


1994 


* 



« Parameters for BASES 



(1) Dimensions of integration etc. 



of dimensions : 
of Wilds : 

of sample points 
of subregions 
of regions 
of Hypercubes 



Ndim = 

Nwild 

Ncall 

Ng 

Nregion = 
Ncube 



5 ( 50 at max . ) 

5 ( 15 at max.) 
9375(real) lOOOO(given) 
50 / variable 

5 / variable 
3125 



(2) About the integration variables 











+- 




-+- 




i 


XL(i) 




XU(i) 


+- 


IG(i) 




Wild 


1 


.OOOOOOE+00 


1 


OOOOOOE+00 




1 


-+- 


yes 


2 


.OOOOOOE+00 


1 


OOOOOOE+00 




1 




yes 


3 


.OOOOOOE+00 


1 


OOOOOOE+00 




1 




yes 


4 


.OOOOOOE+00 


1 


OOOOOOE+00 




1 




yes 


5 


.OOOOOOE+00 


1 


OOOOOOE+00 




1 




yes 










-+- 




-+- 





(3) Parameters for the grid optimization step 
Max.# of iterations: ITMXl = 7 
Expected accuracy : Accl = . 1000 '/. 



(4) Parameters for the integration step 



Max.# of iterations: 
Expected accuracy : 



ITMX2 
Acc2 



15 
.0500 '/. 



« Computing Time Information 

(1) For BASES 
Overhead 

Grid Optim. Step 
Integration Step 
Go time for all 



H: M: Sec 
0: 0: 7.46 
0:13:41.90 
0:29:13.41 
0:43: 2.77 



(2) Expected event generation time 
Expected time for 1000 events : 



17.81 Sec 
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Date: 97/ 7/ 1 11:49 

Convergency Behavior for the Grid Optimization Step 



<- Result of each iteration -> <- Cumulative Result -> < CPU time > 
IT Eff R_Neg Estimate Acc '/. Estimate (+- Error ) order Acc "/, ( H: M: Sec ) 



1 


100 


.00 


1 


351E- 


-01 


5 


295 


1 


351125(-^- 


.071544)E-01 


5.295 





1 


57 


29 


2 


100 


.00 


1 


456E- 


-01 


1 


470 


1 


447336 (-^- 


.020502)E-01 


1.417 





3 


54 


83 


3 


100 


.00 


1 


430E- 


-01 


1 


211 


1 


437373(-^- 


.013232)E-01 


.921 





5 


52 


58 


4 


100 


.00 


1 


448E- 


-01 


1 


130 


1 


441760 (-^- 


.010290)E-01 


.714 





7 


50 


20 


5 


100 


.00 


1 


453E- 


-01 


1 


205 


1 


444741(-^- 


. 008872) E-01 


.614 





9 


47 


20 


6 


100 


.00 


1 


457E- 


-01 


1 


171 


1 


447447 (-^- 


.007871)E-01 


.544 





11 


44 


74 


7 


100 


.00 


1 


452E- 


-01 


1 


264 


1 


448097(+- 


.007233) E-01 


.499 





13 


41 


90 



Date: 97/ 7/ 1 11:49 
Convergency Behavior for the Integration Step 



<- Result of each iteration -> <- Cumulative Result -> < CPU time > 
IT Eff R_Neg Estimate Acc Estimate (-f- Error ) order Acc '/. ( H: M: Sec ) 



1 


100 


.00 


1 


467E- 


-01 


1 


159 


1 


467078 (-^- 


.017002)E- 


-01 


1.159 





15 


39 


79 


2 


100 


.00 


1 


436E- 


-01 


1 


111 


1 


450578(-^- 


.011636)E- 


-01 


.802 





17 


37 


65 


3 


100 


.00 


1 


469E- 


-01 


1 


101 


1 


456866 (-^- 


. 009447) E- 


-01 


.648 





19 


35 


01 


4 


100 


.00 


1 


450E- 


-01 


1 


119 


1 


455106(-^- 


. 008164) E- 


-01 


.561 





21 


31 


73 


5 


100 


.00 


1 


453E- 


-01 


1 


100 


1 


454677 (-l- 


.007271)E- 


-01 


.500 





23 


28 


23 


6 


100 


.00 


1 


457E- 


-01 


1 


158 


1 


454980(-^- 


. 006677) E- 


-01 


.459 





25 


24 


70 


7 


100 


.00 


1 


450E- 


-01 


1 


146 


1 


454232 (-^- 


. 006195) E- 


-01 


.426 





27 


21 


81 


8 


100 


.00 


1 


446E- 


■01 


1 


091 


1 


453082(■^- 


.005766)E- 


-01 


.397 





29 


18 


61 


9 


100 


.00 


1 


484E- 


■01 


1 


199 


1 


455997 (■^- 


.005485)E- 


■01 


.377 





31 


15 


24 


10 


100 


.00 


1 


443E- 


■01 


1 


202 


1 


454854(■^- 


.005230)E- 


■01 


.359 





33 


12 


25 


11 


100 


.00 


1 


453E- 


■01 


1 


195 


1 


454674 (■^- 


.005008)E- 


■01 


.344 





35 


8 


71 


12 


100 


.00 


1 


481E- 


■01 


1 


212 


1 


456559(■^- 


. 004824) E- 


■01 


.331 





37 


5 


18 


13 


100 


.00 


1 


457E- 


■01 


1 


248 


1 


456616 (■^- 


.004662)E- 


■01 


.320 





39 


2 


19 


14 


100 


.00 


1 


462E- 


■01 


1 


212 


1 


456987(■^- 


.004509)E- 


■01 


.309 





40 


58 


71 


15 


100 


.00 


1 


464E- 


■01 


1 


360 


1 


457306 (■^- 


.004398)E- 


■01 


.302 





42 


55 


31 
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Histogram (ID = 1 ) for 
X d(Sigma)/dx 
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4 


.491+- 
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3 
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I 


4, 


.467+- 
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E 
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E 
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4, 
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4, 


.453+- 
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E 
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4, 
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I 


4, 
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E 
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4, 


.500 


I 


4, 


.469+- 


.080 


E 


I 


4, 


.650 


I 


4, 


.402+- 


.077 


E 


I 


4, 


.800 


I 


4, 


.381+- 


.082 


E 


I 


4, 


.950 


I 


4, 


.554+- 


.084 


E 


I 


5 
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I 


4 
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.076 


E 
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Number of generated events 

Generation efficiency = 

Computing time for generation = 

for Overhead = 

for Others = 

GO time for event generation = 
Max. number of trials MXTRY 

Number of miss-generation = 



10000 
23.794 Percent 

443.281 Seconds 
8.240 Seconds 
10.349 Seconds 
461.870 Seconds 

50 per event 

51 times 
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Original Histogram (ID = 1 ) for Momentum of Particle 3 
Total = 10000 events "*" : Drig. Dist. in Log Scale. 

X d(Sig/dx) dN/dx l.OE-05 l.OE-04 l.OE-03 
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X d(Sig/dx) dN/dx "0" : Generated Events. ( Arbitrary unit in Log ) 
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